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Abstract 

in 
o 
o 

The influence of continuous measurements of energy with a finite accuracy 
is studied in various quantum systems through a restriction of the Feyn- 
man path-integrals around the measurement result. The method, which is 
equivalent to consider an effective Schrodinger equation with a non-Hermitian 
Hamiltonian, allows to study the dynamics of the wavefunction collapse. A 
numerical algorithm for solving the effective Schrodinger equation is devel- 
oped and checked in the case of a harmonic oscillator. The situations, of 
Q physical interest, of a two-level system and of a metastable quantum- well 

O are then discussed. In the first case the Zeno inhibition observed in quan- 

tum optics experiments is recovered and extended to nonresonant transitions, 
in the second one we propose to observe inhibition of spontaneous decay in 
mesoscopic heterostructures. In all the considered examples the effect of the 
continuous measurement of energy is a freezing of the evolution of the system 
proportional to the accuracy of the measurement itself. 
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I. INTRODUCTION 



Technological progress in the measurement of physical quantities, expecially in quantum 
optics, in the physics of superconducting coherent devices such as SQUIDs, in the study 
of mesoscopic structures and in experimental gravitation, has revived the fifty-years old 
debates on foundations of quantum theory [1,2]. Collapse of the wavefunction, repeated 
measurements on a single object, inhibition of transitions under the effect of measurements 
are now also the language used by the experimentalists involved in this so-called quantum 
phenomenology. Ideal experiments which were thought in the thirties to sharp the para- 
doxical aspects of quantum theory when confronted to classical physics are now feasible. 
This originates a demand for quantitative predictions of quantum measurement theory. In 
general, among the observables of interest for the experimenters a crucial role is played by 
canonical coordinates and energy. The first kind of observable has been discussed exten- 
sively until now, expecially in connection to the fundamental limitations in the accuracy 
of position measurements [2,3]. In this paper we discuss a quantum measurement model 
which allows to describe the general features of an energy measurement independently of 
the particular measuring apparatus used to perform it. The model is based upon restriction 
of path-integrals around the measurement result. After a general description of the tech- 
nique in Section II we analyze the case of a harmonic oscillator in Section III. This allows 
to check the numerical technique on a well established potential. The analysis turns out to 
more physical situations in Sections IV and V. In Section IV the case of a two-level system 
is analyzed and its link to the quantum Zeno effect is estabilished. In Section V the case 
of a properly designed heterostructure is analyzed to study the feasibility of an experiment 
exploiting spontaneous decay of localized electron states. This allows also for a comparison 
between the path-integral model of energy measurement and a more concrete one taking 
into account the detailed properties of a meter. In Section VI the conclusions are given: 
inhibition of the evolution of a system is obtained whenever measurements of energy are 
performed on it in the quantum regime of sensitivity. 



II. GENERAL FORMALISM 

As a consequence of the interaction between the measured and the measuring systems, 
decoherentization arises. The phenomenon can be derived in the framework of standard 
quantum mechanics by summing over the degrees of freedom of the measuring system [4]. 
The measured system then evolves influenced by the interaction with the measuring system 
in such a way that the off-diagonal terms of its density matrix get damped. Of course, 
the exact evolution equation for the density matrix of the measured system depends on the 
speciflc nature of the measuring system, which is an inconvenient if general considerations 
are of interest. The path integral formalism with constraints added to force the system to 
move in a space region [5] or along an individual Feynman path [6] is a simple and effective 
tool for describing the essential features of a continuous measurement without reference to 
meter-dependent details. Here we will make use of the restricted path-integral approach 
proposed by Mensky [7,8]. 

Two different situations can arise: a prediction on the evolution of the measured system is 
requested before the continuous measurement is performed; a complete quantum mechanical 
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description of the measured system, namely its wavefunction, is desired after the continuous 
measurement has been performed. In the first case, the a priori analysis^ the outcome of 
the measurement is unknown and a sum over all possible results should be considered. The 
evolution of the measured system (during the measurement) is described by the density 
matrix. In the second case, the a posteriori analysis^ the outcome of the measurement is 
known and a sum over all Feynman-paths compatible with the result should be considered. 
The measured-system evolution is described by the wavefunction. In both cases, the infiuence 
of the measurement on the evolution of the measured system is represented by a restriction 
of the possible choices (density matrix elements or paths) according to a weight functional. 
The strength of the restriction, measuring the amount of disturb per unit time fed by 
the measuring apparatus into the measured system, completely defines the measurement 
operation. 

The a posteriori analysis has been shown to be very effective in understanding the ac- 
curacy of measurements of position in non-linear systems, monitored in a continuous [9] 
and impulsive way [10]. Here we will show how it applies to the case of measurements of 
energy. Let us suppose that a continuous measurement of energy with result E[t) has been 
performed. The probability amplitude for the measured system going from point q at time 
to point q' at time t in the configuration space during the measurement is 

K[E]iq\t;q,0) = Jb]e* /o g]. (1) 

Jg 

The weight functional W[e] restricts the integration around the measurement result. From a 
formal point of view a very simple and useful choice for W[e] is a Gaussian functional 

W[E][q,p] = exp (^-K j\H{t') - Eit'yfdt'^ . (2) 

The constant k has dimensions l/(time X energy^) and can be put in the form = re^. 
The constants r and e have a well defined meaning in relation to the specific nature of the 
measuring apparatus. Konetchnyi, Mensky, and Namiot [11] have demonstrated that in the 
case of a continuous measurement of position the corresponding weight functional 

W[4q] = exp (^--^ J\q{t') - x{t')fdt'^ (3) 

coincides with that obtained for a measuring apparatus made by an infinite set of quantum 
oscillators, each interacting with the measured system within the length scale A and for a time 
interval r. A similar calculation in the case of a continuous measurement of energy allows us 
to interpret r and e as the characteristic time and energy scales of the interaction between 
the measuring and the measured systems. In other words r and e are the duration and the 
error of each microscopic measuring event. An interpretation in terms of the macroscopic 
measurement operation is also possible. If the measurement lasts a time At = Nt we 
can write = re^ = AtAE^ where AE = ej^/N. Here AE represents the error of 
the continuous measurement and it decreases, as expected, with the squared root of the 
measurement time. 

Due to the choice of the weight functional, the probability amplitude can be written in 
the form 
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'q 

where we introduced the effective Hamiitonian 



H.ff = H-—(H-Er. (5) 



re2 



Due to the presence of a fourth power of the momentum in Hej / the functional integral (4) 
is meaningless from a rigorous point of view and we will use it only in a formal way to get 
a differential equation of motion. During the measurement the evolution of the system is 
given by the modified Schrodinger equation 

d 

ih—'ijj{x,t) = Heff'ip{x,t). (6) 

We note that the effective Hamiitonian (5) is not Hermitian. However, the eigenstates of 
the Hamiitonian H of the unmeasured-system are also eigenstates of (5). We have just a 
modification of the time evolution of these eigenstates with a damping proportional to the 
difference between the corresponding eigenvalue and the measurement result. This allows 
us to dynamically describe the process of an incomplete and/or not instantaneous collapse 
of the state of the measured system. 

Let us consider the case of H with a discrete energy spectrum 

Hcf>.,{x) = E.,cf>4x). (7) 

Let ip[x^O) be the state of the system at the beginning of the measurement. We expand this 
state in the base 

V'(x,o) = ^c4o)</)4x) (8) 

n 

with the normalization J2n — 1- By inserting (8) in (6) we have an evolution equation 

for the coefficients c„(t) 



d i 



En - ^{En - E{t)f 



Cnit) (9) 



with solution 



At the generic time t during the measurement the measured system is completely described 
by the wavefunction 

ij{x,t) = J2cn{t)M^)- (11) 

n 

Due to the presence of the anti-Hermitian part of Hejj^ ip{x^t) loses its initial normalization 
according to the selection rule imposed by the measurement result. The probability to have 
the measured system in the eigenstate n at time t is 
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At the end of the measurement t = At, the probabihties P„ have a simple interpretation 
in terms of the total error AE of the measurement. When AE ^ \En — E\ (for simplicity 
we consider the case of a constant measurement result E in the following discussion) the 
amount of back-action of the meter on the system is negligible and the probability to have 
the system in the eigenstate n is close to that in absence of measurement. However, the 
above disequality can not be valid for all n and a generic state of the system, i.e. containing 
contributions from all n, is always disturbed by the measurement at some extent. If the 
measurement result is some definite eigenvalue E = Ei, the wavefunction of the measured 
system approaches the eigenstate (f)i as much closely as much AE gets smaller. In the limit 
of an infinitely precise measurement a complete collapse of the wavefunction is obtained. 



III. MEASUREMENTS IN A HARMONIC OSCILLATOR 

Here we apply the general formalism introduced in the previous section to the case of a 
harmonic oscillator described by the Hamiltonian 

This simple system allows us to show with analytical tools the mechanism of the continuous 
collapse due to an energy measurement. It allows us also to test the numerical method 
described in appendix A that we will use in the analysis of more complicated systems. 

The evolution of a wavefunction collapsing during an energy measurement is shown in 
Fig. 1. We consider an initial Gaussian state 

1 x2 

^{x,{)) = ^^=e~^ (14) 



which allows a simple analytical evaluation of ip[x^t) in terms of the harmonic oscillator 
eigenstates. The analytical results differ for less then 0.1 % from the numerical results 
shown in Fig. 1 obtained with the algorithm of Appendix A. In absence of measurement, 
i.e. AE = oo, the wavefunction remains a Gaussian with its width oscillating with period 
tt/lo. When the measurement is effective only the projection of ip[x^t) on the eigenstate 
with energy closest to the measurement result survives. The collapse is gradual and only 
for very precise and/or long lasting measurements it is fully accomplished. This behavior 
constitutes the basis of an effective relaxation method for computing numerically eigenvalues 
and eigenstates of Schrodinger operators [12]. 

Beside its simplicity, the harmonic oscillator is quite inappropriate for describing the 
interplay between energy measurements and stimulated transitions, a situation of wide in- 
terest in the applications, expecially in quantum optics. Indeed, due to the constant spacing 
of the spectrum, under a forcing term at resonance all the levels are occupied and no sta- 
tionary regime for the stimulated transitions among the levels can be obtained. Systems 
with nonuniform level spacing or, for simplicity, two-level systems should be considered. 
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IV. STIMULATED TRANSITIONS IN A TWO-LEVEL SYSTEM 



The interplay between the transitions stimulated in a system by an external perturbation 
and the effect of a continuous measurement of energy have been recently discussed within 
the present formalism [13]. It turns out that the so-called quantum Zeno paradox introduced 
in [14] and [15,16] and observed in [17] is a particular example of the influence of the meter 
on the measured system (see also [18] for the debates following [17]). Here we discuss the 
differences between on-resonance and off-resonance perturbation. The case of off-resonance 
transitions is important expecially for schemes of measurements of energy in electromagnetic 
cavities based upon dispersive techniques recently proposed [19]. 

Let us return to the general formalism of Section II. Transitions among the levels of H 
are obtained under the action of an appropriate external perturbation V{t) which is added 
to the effective Hamiltonian of Eq. (5). The decomposition of the state in terms of the 
eigenstates of the unmeasured and unperturbed system can be used again. The evolution 
equation for the coefficient c„ contains also a term proportional to the perturbation strength 
and nondiagonal in the index n, 



d i 



ih 2 



En - , , A En - Eit)) 



^it) - ^Y.^nm{t)c^{t) (15) 



where Km(0 = 1^(0 1*/*™)- 

A particularly simple picture is obtained for a two-level system with energies Ei and £"2- 
Assuming a perturbation potential Vii = V22 = and V12 = V21 = Voe"^^*"*'') with Vq real, 
the solution of the system (15) is 



Ei^ {Er-E) 



2 



Ci(t) = exp <^ ~i~t A j_ A ^ "I" 



h AtAE 



gci(O) + e''-*>C2(0) . 



Ci(0) cos [wt) H : sin {wt] 

iw 



(16) 



\ E2 

C2[t) = exp < 



:e2 - Ef 

AtAE^ 



t — iqt 



€2(0) cos (wt) 



gc2(0) 



-ituto 



pciiO) 



sin [wt] 



iw 



[17) 



where p = Vo/h, q = SE/2h + with = [{E2 - Ef - [Ei - Ef]l2AtAE^ and w = 
\/ -|- p^. The resonance condition is measured by the parameter 

5E = huj- {E2 - El). (18) 

In order to evidence the Zeno effect in a specific example let us suppose that initially 
the system is in the state 1 and the result of the continuous measurement is = £"1. The 
probability Pi{t) to have the system at time t < At in the state 1 is 



\ciitW 



1 



|ci(t)|2+|c2(t)|2 1^ 



Vo 



n.Q-i5E/2+n.w cot(wt) 



(19) 



where 
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hn - -se] 



2 J 



(20) 



and 







(21) 



2AtA^2 



When the measurement error is large, i.e. — > 0, the system oscillates between levels 1 
and 2 with Rabi frequency \J^yQ + SE'^ /h. Complete transitions to level 2, i.e. Pi = 0, 
are obtained only with a resonant perturbation 5E = 0. In the opposite limit of accurate 
measurements, when w is imaginary, an overdamped regime is achieved in which transitions 
are inhibited. A critical damping is observed when w is minimal, i.e. at a measurement 
error 



which defines the borderline between the Rabi-like behaviour [AE > AEcrit) and the Zeno- 
like inhibition [AE < AEcrit) in terms of the instrumental accuracy of the meter. The 
behavior of Pi{t) is shown in Fig. 2 for 5E = and in Fig. 3 for 5E ^ for different 
choices of the measurement error AE. In panels a) we show the Rabi-like behavior (solid 
line) in comparison with the corresponding results with no measurement performed (dashed 
line). Fven in this regime the measurement has the effect to slightly decrease the transition 
frequency. In panels c) we show the Zeno-like behaviour (solid line) in comparison with the 
full inhibition occurring for AE — > (dashed line). In panels b) we show the behavior in the 
critical regime just above (solid line) and below (dot-dashed line) AE = AEcrit- Notice that 
in the off-resonance case the oscillations in the Rabi regime do not reach the zero even in 
absence of measurement. Moreover they get damped in presence of measurement even if the 
system is far from the critical region due to the imaginary term in Fq. (20). On the other 
hand no significant difference is observed between on-resonance and off-resonance cases in 
the Zeno regime. There the effect of the measurement is so strong to completely dominate 
any external perturbation. 



The experiment designed to study quantum Zeno effect in atomic spectroscopy has been 
considered by some authors as not representative of the whole paradox because stimulated 
transition is used to connect two different levels. In this section we discuss a proposal for 
observation of quantum Zeno effect in systems subjected to spontaneous emission. The 
system we discuss is schematically depicted in the inset of Fig. 4. Metastable states are 
obtained in a semiconductor heterostructure where a barrier of width L separates a well of 
width / from a collector region. Let / be the point separating the collector region [/, -|-oo[ 
from the well region ] — oo, /]. For simplicity we suppose equal valence-band offsets Vq and 
we neglect electron-electron interaction. Flectrons generated inside the well, e.g. by a laser 
pulse, relax almost instantaneously in the el state (the low lying state of the quantum well) 
and then start tunneling outside the well. The validity of this picture implies a tunneling 





V. SPONTANEOUS DECAY FROM A QUANTUM WELL 
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time much greater than the phonon relaxation time, i.e. a barrier width L not too smaU. 
On the other hand we suppose the tunnehng process faster than radiative and nonradiative 
electron-hole recombinations so that electrons leave the el of the well only via tunneling into 
the collector. In absence of measurements an exponential decay of the charge trapped in 
the well is obtained. As explained in Appendix B, the single- well metastable potential has a 
low lying resonance state (j)i{x) with complex eigenvalue Ei — ^^l/2 which almost coincides 
in the well region with the el state (j)^[x). Electrons are initially in the el state, i.e. their 
wavefunction is ip{x^ 0) ~ 4>i{x) for x < f. The probability to have an electron in the region 
] — oo, /] at time t, the non-decaying probability, is 

P{t)= |</)^(x)e-*(^i-''ri/2)t|2^^ ^ g-ritA_ ^23) 

J — oo J — oo 

Notice that P{t) is related to the electronic charge Q{t) measured in the collector region by 
Q{t) = Ne[l — -P(t)], N being the number of electrons photogenerated in the well and e the 
electric charge. 

Let us suppose that an energy measurement with known result is performed on the 
collected charge. A practical realization of such a kind of meter is obtained by growing a 
second barrier of width separated by a well of width from the measured system (see the 
inset of Fig. 5). The meter system acts as a filter allowing only the charge in a certain energy 
window to enter the collector region. Center and width of the energy window are specified 
by the well width and the barrier width of the meter, respectively. The energy 
measurement modifies the exponential decay of Eq. (23) depending on the measurement 
result and accuracy. We will report a comparison between the description in terms of the 
restricted path-integral method and the description in terms of the above-mentioned specific 
meter. 

Since we are interested to follow in time the probability P{t) (or the collected charge 
Q{t)) without specifying the duration of the measurement, it is convenient to express the 
measurement feedback through the microscopic quantities r and e. Notice that there is only 
one degree of freedom in choosing these constants, the combination re^ only being relevant. 
We put T = h/e leaving e as a free parameter. By assuming that the result of the energy 
measurement is constant E = Ei, the non-decaying probability is 

P(t)~ f |</>-(x)e-^(^i-ri/2)tg-^(i^,-.Ti/2-i^)^t|2^^^g-r.„tA (24) 

J—OO 

where 

Te// = Ti (l - ^) . (25) 

The decay of the electron measured to be inside the well in the el state is frozen propor- 
tionally with the measurement precision. Equation (25) is rigorously valid only for e ^ Fi. 
For e < Fi/2 it gives an absurd negative decay rate, however, already for e ^ Fi the re- 
sults obtained from (25) are in good agreement with the numerical simulations reported in 
Fig. 4. In the numerical simulations P{t) is obtained from the wavefunction evaluated in 
a space-time lattice with the algorithm described in Appendix A. We have chosen / = 60 
A, iv = 50A,Vo = 0.1 eV and an electron effective mass m = 0.1 mo, mo being the bare 
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electron mass. With these parameters the first resonance level evaluated as explained in 
Appendix B is = 3.72 IQ-^ eV and Ti = 6.99 IQ-^ eV. 

Now we turn to the measurement performed with the model meter previously discussed. 
Since the measurement result has been chosen to be = £"1, we have = /. We put also 
Ljn = L. The measured system couples quantum mechanically with the meter system in a 
well known manner. The previous first resonance splits into two resonant states (j)i{x) and 
</)2(x) for which we evaluate Ei = 3.45 IQ-^ eV, Ti = 2.86 10"^ eV and E2 = 4.01 IQ-^ 
eV, r2 = 4.35 10"'* eV. Assuming that the electrons are initially in the el state of the 
unmeasured system, i.e. their wavefunction is ?/'(a;,0) ~ [(j)i{x) + (/)2(a;)]/\/2, for x < /, we 
have 

P{t) c^-f \4>r{x)e-^n(E^-^^^l2)i ^ ^oo(^^g-^(i^.-.T./2)t|2^^ ^ i (-T^i/n ^ g-r.tA\ _ ^26) 

^ J — GO 2 ^ 

Notice that the non-decaying probability includes integration over all the space but the 
collector region, i.e. ] — 00, /] includes the meter. The rough prediction of Eq. (26) is shown 
in Fig. 5 (dashed line) in comparison with the corresponding exact numerical simulation 
(solid line). The exact probability P{t) decays oscillating around the curve of Eq. (26) with 
frequency [E2 — Ei)/h. The model meter considered with barrier width = L is expected 
to represent an energy measurement of precision e ~ Fi ~ r2 and duration r ~ h/Ti. 
Indeed, a comparison of Fig. 4 and Fig. 5 shows that this prediction is confirmed. 

In the example we have discussed the realization of the back-action of the meter, i.e. 
its infiuence on the measured system, is quite evident. As the barrier acting as a meter 
measures the energy of the transmitted electrons with increasing accuracy, the resonances 
get narrower and the refiected current increases, coming back to the measured system. In 
the opposite limit, when the measurement is not accurate, the refiected current and the 
consequent back-action on the measured system are negligible. 

The results obtained within the restricted path-integral formalism both in this section 
and in the previous one depend on the Gaussian ansatz (2) for the weight functional. This 
choice is not only mathematically convenient but corresponds to a specific, still quite general, 
model of measuring apparatus obeying the quantum-dynamical semigroup property [11,20]. 
Therefore qualitative agreement with the results obtained in a real measurement must be 
expected. On the other hand, detailed behaviors in a real measurement could be taken into 
account only by an accurate modellization of the employed meter. A comparison of Fig. 
4 and Fig. 5 helps to clarify this point in a specific example. The restricted path-integral 
prediction correctly reproduces the main feature of the barrier-modelled meter, namely the 
variation of the decay time constant with respect to the unmeasured system case. On the 
other hand, the oscillations of P{t) in Fig. 5 superimposed to the smooth decay curve 
(dashed line) have a period hir / [E2 — Ei) and are the coherent oscillations of the charge 
between the well of the system and the well formed by the barrier meter. Of course these 
oscillations are not obtained in the restricted path-integral approach. This is a desirable 
feature when looking at a general formulation of quantum theory of measurement. 

VI. CONCLUSIONS 

The time evolution of various systems during a continuous measurement of energy has 
been studied within a very general formalism. It has been shown that, regardless of the spe- 
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cific realization of the meter, the measurement causes a partial freezing of the time evolution 
of the system. This fact is already known in quantum optics where a Zeno-like experiment 
has been performed. The generality of our approach allows not only to recover the quantum 
optics case but also to describe quantitatively a proposal for analogous experiments using 
man-made mesoscopic structures. In this case one could have the possibility to test Zeno 
inhibition on systems manifesting spontaneous decay with easier control of the parameters 
of the structure under measurement. 

Our model can be considered as a quantitative dynamical approach to the energy-time 
uncertainty relationship already debated in the literature [21,22]. Let us define the infiuence 
time 6 = h/Tejf — h/Ti expressing the modification of the spontaneous decay time under 
the infiuence of an energy measurement with error e. By using Eq. (25) in the appropriate 
limit of validity e ^ Fi we have 
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h 

2- (-27) 



High precision measurements of energy lead to changes of the evolution time-scales for the 
measured system. In the limit of an ideal, infinite accuracy, measurement the evolution of 
the system is completely frozen and the so called quantum Zeno effect is recovered. 



APPENDIX A: NUMERICAL SIMULATIONS 

Here we briefly discuss the algorithm used for the numerical simulations of Eq. (6). The 
partial differential equation is transformed in a set of algebraic equations by discretizing the 
space-time x — t in an appropriate two dimensional lattice. Let us suppose that we want to 
follow the temporal evolution of the initial state ip{x^ 0) during the time interval [0, tmax]- Let 
[XfTiim Xjnax] be the space interval where the wavefunction ip[x^t) can be considered localized 
for < t < tjnax- The space-time lattice is deflned by putting 

X ^ X.rara+ J^X J = 0, 1 , 2, . . . , J + 1 (Al) 



t^nAt n = 0,1,2,..., iV (A2) 

where the total number of points in the lattice is related to the lattice constants Ax and At 
by 

{/J ~\~ 1)^^ ^raax '^min l'^^) 



NAt = tr^ax. (A4) 
The wavefunction reduces on the lattice to a matrix 

(A5) 

where rp'^ are known. The effective Hamiltonian operator Hejj is 
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H, 



H-—AH-Ef=(l + —E\H-—H'-—E' 



;2 /o„„\a2 /a^2 



(A6) 



and includes the action of the Hamiltonian H = —{h /2m)d / dx + V[x^t) and its squared 
on tlj{x,t). The action of H on tlj{x,t) gives 



Htpix^t) 



(A7) 



where = For we have 



2mAi 



J + 2 



2m Ax"^ 



i^^iV'j+i - 2y/v; + ^-iV'^-i + ^/(V'j+i - 2V'; + V'j-i)] • (A8) 



The discretization of the operator ihd/dt is obtained with the Cayley approximation for 
the exponential time evolution operator 



-j-KAt 



lAt ■ 
1 + -^K 
2h . 



1 



lAt 
~2h 



-K 



(A9) 



which preserves unitarity for K Hermitian and is correct 0{At^) for K time-independent. 
In our case Hpj j depends on time and is the sum of a Hermitian term and an anti-Hermitian 
one. The Cayley approximation is still the right discretization scheme since it avoids the 
stability problems arising from the Hermitian part. The finite difference equation for the 
time evolution is the Crank-Nicholson scheme 



V'(t + At) 



1 



lAt 
~2h 



Heffit) 



(AlO) 



which corresponds in the lattice to the following linear system 



(Aii; 



where 



7 



/5- 



lAt 



2re^ \2mAx'^ ^ 
ihAt 

-47 



67 + 



2mA 
ihAt 
2mAx'^ 



x^ 



l+^ '3 
lAt 



+ 



2h 



V- - —JE - v: 



n\2 



(A12) 
(A13) 
(A14) 



By using the condition V'"i = V'o = V'j+i = V' 
wavefunction is localized in \x 



V+2 



which reflects the fact that the 



mm •} '^max] •} 



we can write Eq. (All) in matrix notation 
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(% + ni+') = (% - i^- ^, J = 1, . . . , J (A15) 
where the matrix 7^" is 5-diagonal 



/«1 


/3r 
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«}-2 


/5j-2 


7 
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■ 7 


/5j-2 








V 
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7 




a} 


J 



(A16) 



Starting from rp'^ the solution of the above system at each time step gives the unknown 
V'""'"^ in terms of the known Notice that the matrix 7^" does not depend on time if the 
potential V is time independent. 

Due to its 5-diagonal nature the linear system in Eq. (A15) is solved by LU decomposition 
accomplished through standard-library subroutines. The number of iterations required for 
the solution grows linearly with the matrix dimension J. The computer time needed for the 
simulation of the full evolution is about C N J , with C ~ 7.6 /is in a VAX 7000-610 machine. 



APPENDIX B: RESONANCE STATES OF THE METASTABLE POTENTIAL 

We want to derive a semiclassical formula for the resonances A„ = — ^^„/2 of the 
metastable square-well potential 



V{x) = Vo l[-oc,a] + Vo l[b,c] + V'o l[d,f] 



iBl] 



where a<b<c<d<f and l[a,b] means 1 for a < x < b and elsewhere. The 
eigenfunctions of the stationary Schrodinger equation 



2m dx"^ 



+ V{x) 



iB2) 



are such that (pnix) oc exp (i\/Kix^ , for x — > oo. We also observe that, as a special case of the 

method of complex scaling [23], </>n(a;) is of class on the contour 7 =] — 00, /] U e'^ [/, -|-cx3[ 
if ^ > is conveniently chosen. In the limit of f — d large, we approximate En with the 
eigenvalues of the confining potential V^{x) = Vq l[-oo,a] + M3 l[6,c] + M3 l[(i,+oo[- The functions 



9n 



+ C-e- 
+ D-e 



ikx 
■kox 
-ikx 



X < a 
a < X < b 
b < X < c 
c < X < d 



(B3) 
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where 



/2m 



(B4) 



and 



2m 

^0 = \/-^(^o - E„ 



are the bound states for the potential V^{x) when F'^ 
points a, 6, c and d gives: 



B-- 



~2ik 



(B5) 



0. The requirement at the 



(B6) 



2ko 



B^e'^\ko ± ik) + B-e-'^\ko T 



(B7) 



= — — \C+e''°'(ik ± ko) + C-e-''°'(ik T ko 
2ik L 



fB8) 



2ko 



D^e'^\ko ± ik) + D-e-'''\ko T ik) 



(B9) 



The eigenvalues are found by solving numerically the equation F~^[Fn) = 0. The nor- 
malization of (j)^[x) is achieved by choosing 



A- 



^2koa 



2ko 



+ {\B+\^ + \B-\^) (6-a) + 23f?e{5+5-*e'''=("+^)} 



+2Ue\C+C-*\{c-b) + \C+\^ 
+ (|i)+|2 + |/)-|2j {d-c) + 2Ue\D+D-*e''^'+'^} 



*^^k(c+dU sm[k{d-c)] 

k 



k 

-2kob _ ^2koc 



2ka 



-2kod 



2ko 



(BIO) 



Once Fn is known, the corresponding r„ can be found by observing that in the limit of 
f — d large we have the semiclassical approximation 



(l)n{x) = A < 



e^°^ X < a 

C+e'^o^ + C-e-'^"^ b< X <c 
p+^kox ^ p-^-kox d<x < f 



(Bll) 



G+ 



-,ikx 



f <x 



with A, B^, C^, and F evaluated as above and F~^ , obtained through the 
requirement at / 
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F+ = F-^^^ + ^e-^'^o/ (B12) 
ko — ik 

G+ = F--^e-^-f-^^f. (B13) 
ko — ik 

By multiplying the stationary Schrodinger equation (B2) with (/>* and the complex conjugate 
equation with and subtracting the two results, after integration in ] — oo,/"'"], we find 

r„ \U^)?dx = ^2^m \^^^Mn^ ■ (B14) 
j-oo Zm \^ ax J 

By using Eq. (B13) in the r.h.s. and approximating ~ (j)^ in the integral, we get: 

_ Skkl |AF-e-^"/p 

" 2m F + A;2 1 - |AF-e-'=o/|2/2A;o' ^ ' 

The above results apply also to the case of a single-well metastable potential by putting 
a = h = c. 
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FIGURES 

t/At=0 t/At=l/6 t/At=l/3 t/At=l 




FIG. 1. Evolution of the wavefunction collapse under a continuous measurement of energy in a 
harmonic oscillator. Each panel shows, at different times and for different measurement errors, the 
squared modulus of the wavefunction during a measurement with constant result E = E2 = ^hu. 
The initial wavefunction is a Gaussian of width a = and we put At = Stt/cj, h = 2m = 1, 

= 5 • 10"'*. In the right-most panels corresponding to the time t = At the dashed line indicates 
the squared modulus of the eigenfunction (^2 towards which the initial state collapses for a finite 



16 




FIG. 2. Probability Pi{t) to observe a two-level system in the state 1 during a measurement 
of energy and under the influence of a resonant perturbation. The system is in the state 1 at 
time t = and the result of the measurement is E = Ei constant. The measurement error is 
AE = CO (dashed line) and AE = 2AEcrit (solid line) in panel a), AE = 1.07 AEcrit (solid line) 
and AE = 0.99AEcrit (dot-dashed line) in panel b) and AE = O.bAEcrit (solid line) and AE = 
(dashed line) in panel c). We put At = lOirJi/Vo and E2 — Ei = Vq = fi = 1. 
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FIG. 3. As in Fig. 2 but for a nonresonant perturbation with 8E = {E2 — Ei)/2. 
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FIG. 4. Non-decaying probability P{t) for electrons generated inside the metastable well shown 
in the inset during a continuous measurement of energy achieved by restriction of the Feynman 
paths. The measurement result is the energy of the resonant state and we put t = Ji/e for different 
values of e. 
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FIG. 5. Non-decaying probability P{t) (upper solid line) for electrons generated inside the same 
metastable well of Fig. 4 but during a continuous measurement of energy achieved by the model 
meter shown in the inset (dashed line). The lower solid line is obtained in absence of measurement. 
The dashed lines are the results of Eq. (23) (lower line) and Eq. (26) (upper line). 
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